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■ The light-cone wave function of the pion is calculated within the Nambu-Jona- 
Lasinio model. The result is used to derive the pion electromagnetic form factor, 

\ charge radius, structure function, 7r-7~transition form factor and distribution am- 

' plitude. 
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O ' 1 Introduction 

{iLii A light-cone (LC) wave function is a localized (i.e. normalizable) stationary so- 

lution of the LC Schrodinger equation idrl^ij)) = /7lc|^(t)), which describes 
the evolution of a state |\I/(t)) in LC time r = x+ = x° + x^, conjugate to the 
LC Hamiltonian Hi^c = P~ = P^ — P^ [1)2]. In quantum field theory, a station- 
I ary solution has a LC Fock expansion which in principle does not termi- 

te I nate, as can be seen e.g. for the pion, |7r) = ip2\<i<i) + i^sllQg) + i^illQQQ) + ■ ■ ■■ 

The LC wave functions are the amplitudes ipn to find n particles with mo- 
menta Pi (quarks q, antiquarks q or gluons g) in a pion of momentum P, 
i.e. ipnipi, ■ ■ -Pn', P) = (^l^r) = (pi, . . . ,p„|7r(P)). We shall see in a moment 
that the momentum dependence of the ipn is of a very peculiar nature. 

Upon choosing LC rather than ordinary time one benefits from the following 
virtues. First of all, the field theory vacuum becomes 'trivial', in particular 
P~\0) = 0. As a result, disconnected vacuum contributions do not mix with 
the LC wave functions: {0\P'\n) = 0. Unlike in ordinary 'instant-form' quan- 
tization, boosts become kinematical (interaction independent). This is closely 
related to a twodimensional Galilei invariance which leads to a separation of 
center-of-mass and relative coordinates. These features imply that the ipn 
only depend on the frame independent relative momenta Xi = pt/P'^ and 

k±i = XiP± - p_Lj. 
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As the number of difficulties is always conserved, the list of advantages goes 
along with a list of problems. There is the conceptual question how nontrivial 
condensates can arise in a trivial vacuum. The answer is crucial for understand- 
ing spontaneous symmetry breaking in the LC framework. Instead of boosts, 
rotations and parity become dynamical and thus complicated. There are also 
more technical problems. Due to the lack of explicit covariance and rotational 
invariance, renormalization becomes very difficult beyond one loop. There are 
simply not enough manifest symmetry principles providing a guideline. As a 
result, there is an abuncance of counterterms which can even become non- 
local (~ l/^"*")- In addition, LC field theories are constrained systems. This 
makes their quantization nonstandard and somewhat involved. Finally, the 
LC Schrodinger equation in general represents an infinity of coupled, nonlin- 
ear integral equations for the amplitudes ipn- To solve them one has to resort 
to truncations. It should be noted that for an ab-initio calculation of LC wave 
functions basically all this problems have to be solved. 

If it is so hard, then why should one calculate LC wave functions? It is worth 
the effort because their determination implies the knowledge of the entire 
hadron structure (just like the determination of the Coulomb wave functions 
implies the knowledge of the hydrogen structure). Unlike for nonrelativistic 
systems, however, there is a subtlety. The (rclativistic) LC wave functions 
depend on a resolution scale Q, ipn = V'n(Q)- It is well known that there 
are two basic regimes, a 'soft' and a 'hard' one. In the hard regime, one has 
Q Aqcd so that perturbative QCD is applicable. The LC wave functions 
describe the distribution of partons. If the resolution scale Q is of the order 
of Aqcd, in the soft regime, perturbation theory no longer works. Hadrons 
arc believed to consist of cfTcctivc constituent quarks which should again be 
reflected in the nature of the LC wave functions. Low-energy effective held 
theories seem to provide a reasonable description of the soft regime, at least as 
far as spontaneous chiral symmetry breaking (xSB) is concerned. Chiral quark 
models, in particular, with explicit quark degrees of freedom, should yield soft 
LC wave functions describing hadrons as bound states of a few constituent 
quarks. Such a description should be valid for Q ~ 1 GeV. Whether these 
expectations are true will be examined in this contribution. 



2 Formalism 

It turns out that most of the difficulties mentioned in the introduction can be 
circumvented using a particular chiral quark model originally due to Nambu 
and Jona-Lasinio [3]. For two flavors the Lagrangian is 




2 



The model is not renormalizable whence it requires a cutoff A, the value of 
which is fixed by phenomenology. Beyond a critical coupling Gc one has spon- 
taneous xSB together with dynamical mass generation, tuq — > m ~ (0| |0). 
The current quarks thus become constituent quarks with a mass proportional 
to the chiral condensate. As a constituent picture is realized one expects a 
truncation of the LC Schrodinger equation to work reasonably well. We avoid 
to solve complicated constraint equations [4-6] by using a Schwinger-Dyson 
approach as pioneered by 't Hooft [7]. It consists of two steps. First one solves 
the Schwinger-Dyson equation for the quark self-energy S in mean-field ap- 
proximation (equivalent to the large- A^^c limit). This yields the gap equation 
S = const = m (0 1 |0) a- To make sense of the quadratically divergent 
expression for the condensate one uses an invariant-mass cutoff A satisfying 
Mg = {k^ + w?)/x(l — x) < A^, X = k+/A (for an alternative, see [8]). It 
regulates the divergence of the condensate both for x 0, 1 and k± ^ oo, 
and is related to the 3- vector cutoff A3 > |k| by A^ = 4(A3-|-m^). Accordingly, 
the result for the condensate agrees with the usual one [2,3]. 

Knowing the effective quark mass m one performs the second step in the 
program and solves the Bethe-Salpeter equation in the pseudoscalar channel 
(using ladder approximation/large-A^c*). The LC wave function of the pion is 
obtained by three-dimensional reduction (integrating over LC energy). In the 
chiral limit, the result is 



where we have defined k± = kii:ik2 corresponding to the Lz = ±1 components 
of the wave function. A^ is a normalization constant to be determined later. 
The step function implements the invariant-mass-cutoff. It is only for this 
cutoff that the wave function becomes normalizable. 



3 Results 

With the pion wave function at hand one can go on and calculate observables. 
To make life simple I will always work in the chiral limit = = M^. Due to 
this limit the wave function (2) does not depend on x (apart from the cutoff). 
Further simplifications arise in the large-cutoff limit (LCL). This amounts to 
keeping only the leading order in = / h?. As a result one can find analytic 
expressions for all observablesQ. The LC wave function (2) simplifies to '?/'2TT = 

^ The same logic is used in the context of the instanton model [9] where corre- 
sponds to the 'packing fraction' of instantons in the vacuum. 




2mk- ra^ — k\ 
— —2mk^ 



9iA' - Ml) , (2) 
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'ip2ii = and ^b2]l = — i-'2iT — NOi^A"^ — Mq). The nontrivial components of 
the wave function thus just become step functions [10]. We are thus left with 
two parameters, N and A, which are determined by normahzing 1/J2 to unity 
{enforcing a constituent picture) and using the pion decay constant f^^. The 
latter is basically given by the wave function 'at the origin' [11]. This results 
in = a/S/Ztt and A = 47r/7r — 1.16 GeV. The value for A exactly coincides 
with the Gcorgi-Manohar scale [12] below which chiral effective theories make 
sense. It is also consistent with the constraint on the wave function stemming 
from 7r°^27 [11], 

I dx ^p^n 0^) = V?>IU = N . (3) 

With all parameters determined one can use the Drell-Yan formula as dis- 
cussed e.g. in [11] and calculate the pion charge radius. The result, = 
3/47r^/^ ~ (0.60 fm)^, agrees with those obtained in a covariant Bethe- 
Salpeter approach [13] and the instanton model [9]. The slight difference com- 
pared to the experimental value of 0.66 fm stems from the use of the LCL. 
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Fig. 1. left: pion electromagnetic form factor squared vs. momentum transfer 
= q\\ dashed line: pole fit = l/(l+g^r^/6) with r^r = 0.60 fm; full line: same 
formula with experimental value as extracted from the data [14]; right: ttj transition 
form factor times momentum transfer Q^; dashed line: NJL pole formula (4); full 
Une: same formula with replaced by A^/2; full horizontal line: asymptotic value 
2/^; circles: CELLO data [15]; squares: CLEO data [16]. 

A pole fit to the form factor data of Amendolia et al. [14] is shown in the left- 
hand part of Fig. 1. The momentum space wave function (2) readily yields 
the r.m.s. transverse momentum, {kD — A^/10 ~ (370 MeV^), which makes 
it obvious that the pion is a highly relativistic system. 

The valence quark distribution f" (or pion structure function) is given by the 
square of the wave function integrated over k± resulting in f^{x) = 6x{l — x). 
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This coincides with the resuhs of [5,17] and quahtatively agrees with the 
empirical parton distributions of [18] if the relevant scale is defined as = 
{x){l - {x))A^ ~ (600 MeV)2, the mean {x) being 1/2. 

The transition form factor for the process 77* — > 7r° can be calculated accord- 
ing to [11]. One finds the pole formula 

1 1 



which is displayed in the right-hand part of Fig. 1. The low-Q^ behavior is 
fixed by (3) and thus fine. The behavior for — >^ cxo is off by a factor of two, 
which, however, is not bothersome as this regime is way beyond where the 
NJL model makes sense. One can fix the large-Q^ behavior by introducing an 
effective, spin-averaged wave function [10]. This simply amounts to replacing 
N 2N and in (2) and (4). In this case, (3) cannot be maintained 

and the logic based on low-energy chiral dynamics does no longer apply. 

The pion distribution amplitude 0njl is given by the /c^-integral of the wave 
function. As the square of a step function is again a step function the distribu- 
tion amplitude coincides with f" and hence with the asymptotic distribution 
amplitude, 0NjL(a;) = 6x{l — x) = (j)as{x)- Recently this quantity has been 
measured at Fermilab [19]. At a rather low momentum scale of Q — 3 GeV cf) 
turns out to be close to asymptotic. This provides further evidence that the 
approach presented above makes sense. 



4 Conclusions 



I have analytically determined the LC wave function of the pion within a sim- 
ple field theoretic model (NJL) which is known to yield a good description of 
xSB. With the pion wave function at hand I was able to calculate a number of 
observables which were accurate to within 10 %, consistent with the estimated 
limitations of the LCL employed. Particularly important is the result that a 
constituent picture does make sense: higher Fock states seem to be unimpor- 
tant. If one wants to go beyond the LCL, one faces a fine-tuning problem in e^. 
Delicate numerical fits will thus be necessary. It would be interesting to refine 
the presented method by choosing more realistic models. Possibilities would 
be (i) the instanton vacuum [9] which is more closely related to QCD, (ii) an 
elaborate Schwinger-Dyson approach with more realistic kernels and propa- 
gators [20], or (iii) an effective field theory approach in terms of constituent 
quarks which, at the moment, is still in its infancy [12]. 
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